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FOREWORD 


'  This  report  describes  the  research  effort  of  the  Data 
Systems  Division  of  Litton  Industries.  Inc.,  under  Con¬ 
tract  No  pAACi05-70-C-0328.  with  the  L'.S.  Depart- 
Tneni  of  the  Army.  Frankford  Arsenal.  The  objecii'  e 
was  to  perform  a  parametric  study  of  the  Advanced 
Forward  Area  Air  Defense  (AFAADS)  Weapon 
System. 

The  original  RFP.  to  which  this  contractual  effort  is 
responsive,  was  titled  'Parametric  Study,  Advanced 
Forward  Area  Air  Defense  Weapons  System 
(AFAADS),'  and  called  for  studies  specific  to  an  anti¬ 
aircraft  gun  system  based  on  a  weapon  havir.o  the 
same  general  characteristics  as  the  37mm  Gatling  Gun. 

In  the  intervening  time,  since  the  3  May  1968  date 
of  issue  of  the  RFP.  the  Army  generalized  the  use  of 
the  acronym  ‘AFAADS’  to  include  all  elements  of  the 
forward  area  air  defense  of  the  Army,  including  mis¬ 
siles.  At  the  present  time,  it  appears  that  the  effort 
contracted  for  could  more  properly  be  called 
‘LOFADS,’  Low  Altitude  Forward  Area  Air  Defense. 


However,  in  order  to  maintain  consistency  with  the 
wording  of  the  contract,  the  term  AFAADS  has  been 
used  uniformly  in  the  present  report.  The  reader 
should  remain  aware  of  the  fact  that  we  are  dealing 
here  only  with  the  antiaircraft  gun  system,  which  in 
turn  is  an  element  of  the  complete  system  for  the 
forward  air  defense  of  the  Army. 

The  report  is  included  in  two  volumes.  Volume  I 
comprises  the  analyses  and  conclusions  derived  from 
an  in-depth  study  by  Mr.  Herbert  K.  Weiss.  Also 
provided  therein  are  descriptions  of  the  AFAADS 
System  Concept  as  well  as  characteristics  of  AFAADS 
and  its  environment. 

Volume  II  describes  a  computer  simulation  model 
used  to  provide  the  capability  of  predicting  the  per¬ 
formance  of  proposed  (or  existing)  antiaircraft  gun 
systems  against  an  aircraft  in  a  one-to-one  combat 
situation.  Results  of  this  effort  by-Mr.  Martin  P.  Gins* 
berg  were  used  to  derive  recommended  configurations 
and  conclusions  as  presented  in  Volume  I. 
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SECTION  1 
INTRODUCTION 


The  purpose  of  the  AFAADS  computer-simulation 
model  is  to  provide  the  capability  of  predicting  the 
performance  of  proposed  (or  existing)  antiaircraft  gun 
systems  against  an  aircraft  in  a  one-to-one  combat 
situation. 

The  model  is  not  particular  to  one  gun  system  or 
aircraft  but  has  the  ability  to  model  a  variety  of 
systems.  The  modular  design  of  the  simulation  allows 
each  capability  to  be  progressively  refined  as  more 
detailed  definition  of  an  element  of  the  system  is 
developed.  The  simulator  at  present  includes  the  facili¬ 
ties  to: 

a.  Model  realistic  aircraft  paths. 

b.  Model  sensor  behavior  including  the  effects  of : 

(1)  Sensor  lag. 

(2)  Regenerative  tracking. 

(3)  Time-correlated  sensor  noise  and  aim  wander. 

(4)  Data  smoothing. 

c.  Model  the  following  prediction  methods: 

(1)  Linear. 

(2)  Linear  plus  energy  corrections. 

(3)  Quadratic. 

(4)  Quadratic  plus  energy  corrections. 

(5)  Defense  of  a  known  point. 

(6)  Simple  polar. 

d.  Model  the  ballistics  of  any  shell  with  trajectory 
characteristics  suitable  to  representation  as  a 
solution  of  the  two  parameter  differential 
equation: 

dV  =  -kVGdt  (1.1) 

e.  Couple  the  selected  prediction  algorithm  and  the 
selected  ballistics  in  the  time  of  flight  and  shell 
effectiveness  computation?. 

f.  Model  gun  motion  including  the  effects  of: 

( 1 )  Gun  lag. 

(2)  Regenerative  tracking. 

g.  Compute  miss  distance  at  the  time  of  closest 
approach  of  the  shell  to  the  aircraft. 

h.  Generate  periodic  and/or  summary  reports  of 
the  output  from  the: 


(1)  Sensor. 

(2)  Data  smoother. 

(3)  Predictor. 

(4)  Miss-distance  calculation. 

i.  Generate  a  summary  report  including  the  follow¬ 
ing  items  in  any  combination: 

(1)  The  single  shot  kill  probability  of  shots 
equally  spaced  in  time  over  the  duration  of 
the  aircraft  course. 

(2)  The  expected  number  of  killing  hits  for  each 
burst  of  a  series  of  equally  spaced  disjoint 
bursts  over  the  duration  of  the  aircraft  course. 

(3)  The  probability  that  th.  aircraft  is  ‘killed’  for 
each  burst  of  a  series  of  equally  spaced  dis¬ 
joint  bursts  over  the  duration  of  the  aircraft 
course. 

(4)  The  number  of  bursts  of  a  series  of  equally- 
spaced  disjoint  bursts  over  the  duration  of  the 
aircraft  course  which  have  kill  probabilities  in 
excess  of  each  of  five  reference  probabilities. 

(5)  The  number  of  bursts  of  a  series  equally- 
spaced  disjoint  bursts  over  the  duration  of  the 
aircraft  course  which  have  miss  distances  not 
exceeding  each  of  five  reference  distances. 

(6)  The  integral  of  the  single-shot  kill  probability 
over  the  duration  of  the  aircraft  course. 

The  computer  program  is  written  in  basic  FOR¬ 
TRAN  IV  and  is  currently  being  executed  on  the  RAX 
time-sharing  system.  No  timing  estimates  are  available 
as  to  the  actual  running  time,  however,  it  is  less  than 
20  seconds  of  RAX  time  for  a  typical  problem.  Actual 
360/40  central  processor  time  is  probably  much  less 
than  the  RAX  terminal  time. 

The  model  is  structured  in  a  manner  which  main¬ 
tains  the  ‘real  world’  information  pattern  of  the  simu¬ 
lated  gun  system.  This  is  accomplished  by  restricting 
the  access  of  each  simulated  gun  component  to  only 
that  information  which  would  be  available  in  the  ‘real 
world.’ 

The  manner  in  which  information  is  acquired,  gen¬ 
erated.  and  transmitted  through  the  system  is  pre¬ 
sented  in  Figure  1-1.  The  flow  is  as  follows:  from 
aircraft  position  to  the  Data  Acquisition  Event,  which 
fills  the  SENSED  POSITION  File:  to  the  Prediction 
Event,  which  accesses  the  SENSED  POSITION  File 
and  fills  the  PREDICTED  POSITION  File:  to  the  Gun 
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Motion  and  Firing  Event,  which  accesses  the  PRE¬ 
DICTED  POSITION  File  and  fills  the  SHELL  POSI¬ 
TION  File. 

The  preceding  information  flow  pattern  demon¬ 
strates  an  important  feature  of  the  model.  Note  that,  as 
in  the  "real  world.’  any  action  which  the  gun  system 
takes  (or  can  take)  is  ultimately  based  on  information 
acquired  through  its  sensing  system. 

In  ihe  following  sections  the  model  capability  is 
described  in  more  detail.  To  facilitate  the  discussion, 
the  model  is  viewed  in  four  parts;  three  sets  of  func¬ 
tionally  related  program  modules  and  a  communica¬ 
tion  pool  through  which  information  is  transmitted 
between  the  various  modules 

The  first  set  of  modules  is  used  to  control  the  simula¬ 
tion  process.  It  is  composed  of  the  main  program  and 
one  subroutine  which  maintains  a  schedule  of  events 
and  provides  for  central  time  accounting. 

The  second  set  of  modules  is  composed  of  the  Gun 
System  Functional  Event  Routines.  This  is  the  set  of 
routines  which  model  the  various  components  of  the 


gun  system  under  study.  There  are  four  functional 
event  routines: 

a.  Sensor  data  acquisition. 

b  Regenerative  tracking  switch. 

c.  Prediction. 

d.  Gun  motion  and  firing, 

The  third  set  of  modules  is  the  utility  package  These 
modules  are  used  by  the  event  routines  and/or  each 
other  and  may  not  be  invoked  by  the  system  scheduler. 
Report  routines  are  excepted  from  this  rule  since  they 
do  not  change  the  contents  of  any  communication  pool 
variables.  Included  in  this  set  are: 

a.  Aircraft  positioner. 

b.  Time  of  flight. 

c.  Prediction  algorithm. 

d.  Random  number  generator. 

e.  Coordinate  transformation  package. 

f.  Report  generators. 
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Figure  1-1.  Information  Flow  Through  the  AFAADS  Simulation  Model 


SECTION  2 

SIMULATION  CONTROL 


The  simulation  model  represents  the  combat  situa¬ 
tion  under  study  as  a  series  of  events  centered  around 
the  antiaircraft  gun  system.  Event  sequencing  and 
initialization  tasks  are  performed  by  the  master  control 
program  and  the  clock  module.  The  initialization 
phase  consists  of  reading  input  parameters,  perform¬ 
ing  and  storing  the  results  of  once-only  computations, 
and  scheduling  the  first  (a  sensor  data  acquisition 
event)  and  last  (a  report  event)  events  to  be  processed. 
Initialization  tasks  are  accomplished  by  special  mod¬ 
ules  which  are  called  by  the  master  control  program 
prior  to  entering  the  event  sequencing  mode. 

The  logic  of  the  control  process  is  presented  in 
Figure  2-1.  The  particular  order  of  events,  which  will 
occur,  is  determined  by  the  gun  system  and  the  aircraft 
path  under  study.  There  are,  however,  some  principles 
worth  noting:  for  example,  the  first  event  to  occur  will 


always  be  Event  Type  i.  and  it  will  continue  to  occur 
at  its  specified  period  (input  during  initialization 
phase)  until  the  end  simulation  event  With  the  excep¬ 
tion  of  the  end  simulation  event,  no  other  events  will 
occur  unless  they  are  explicitly  scheduled  bv  Event 
Type  1  or  some  other  module,  either  directly  or  indi¬ 
rectly  scheduled  by  an  occurrence  of  Event  Type  1.  The 
occurrence  of  Event  Type  2  indicates  the  end  of  the 
current  sequence  of  simulation  events. 

One  of  two  actions  are  initiated  after  the  occurrence 
of  Event  Type  2.  If  there  are  more  replications  to  be 
performed;  the  event  file  is  cleared,  all  flags  are  turned 
off,  the  initial  Event  Type  I  and  end  simulation  events 
are  scheduled,  and  the  simulation  process  continues.  If 
no  more  replications  are  required,  the  summary  statis¬ 
tics  over  replications  are  printed,  and  the  simulation 
process  is  terminated. 
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SECTION  3 
EVENTS 


3.1  SENSOR  DATA  ACQUISITION  EVENT 

in  the  real  world  all  gun  system  functions,  which 
require  a  knowledge  of  aircraft  position  and/or  rate 
parameters,  acquire  the  information  either  directly  or 
indirectly  from  the  gun  sensor.  In  the  computer  model 
the  sensor  function  is  modeled  by  the  Sensor  Data 
Acquisition  Event.  On  each  entry,  this  event  routine 
stores  current  values  of  sensed  aircraft  position  in  the 
SENSED  DATA  File.  The  gun-system’s  data  smoother 
is  also  modeled  in  this  routine  and  its  outputs, 
smoothed  position  and  rate  values,  are  also  placed  in 
the  data  file. 

In  the  following  discussion,  upper  case  letters  will  te 
used  for  true  values,  primed  upper  case  letters  for 
sensed  values,  underscored  upper  case  letters  for 
smoothed  values  of  position,  and  dot  notation  to  indi¬ 
cate  rates  and  accelerations.  The  time  to  which  values 
pertain  will  be  indicated  by  subscript  expressed  relative 
to  current  time.  Thus,  if  A  is  a  position  coordinate: 

A0  ■*  Current  true  value. 

Aq  m  Current  sensed  value. 

A0  «  Current  smoothed  value. 

AQ  -  Current  smoothed  rate  value. 

A'A(-  Sensed  value  at  current  time  minus  At. 

The  six  dimensional  transformations  from  cartesian 
to  polar,  and  from  polar  to  cartesian  coordinates,  will 
be  represented  by  the  operators;  ‘CTOP’  and  ‘PTOC.’ 
Thus. 

(R,  6, 0,  R,  6,0)  =  CTOP  [X,  Y,  Z,  X,  Y,  Z]  (3.1) 
and 

(X,  Y,  Z,  !)  =  PTOC  |R,0,*,M,0]  (3  2) 

Standard  superscript  arrow  notation  (A)  will  be  used 
to  designate  vectors.  Thus,  the  coordinate  transforma¬ 
tions  become: 

?=  CTOP  [Cj  (3.3) 

and 

r=  PTOC  [Pj .  (3.4) 

On  each  entry  to  the  event  routine,  the  sequence  of 
events  is  as  follows.  The  true  aircraft  position  and  rate 
is  obtained  by  invoking  the  aircraft  position  routine. 
These  values,  which  are  returned  in  cartesian  coordi¬ 
nates,  are  transformed  to  polar  coordinates.  Any  •ad¬ 
justments’  made  for  sensor  lag  and/or  regenerative 
s  *nsor  tracking,  and/or  sensor  noise  are  made  to  the 
polar  coordinates.  The  resulting  sensed  position  in 


polar  coordinates  is  transformed  to  cartesian  coordi¬ 
nates  and  stored  in  the  SENSED  DATA  File.  The  most 
recem  ‘n’  values  stored  in  the  SENSED  DATA  File  are 
used  to  compute  current  values  for  smoothed  position, 
velocity,  and  acceleration.  These  values  are  also  stored 
in  the  SENSED  DATA  File. 

The  compulation  for  sensor  lag  is  based  on  true 
aircraft  rates  since  these  are  the  cause  of  the  phenome¬ 
non  in  the  real  world.  The  lag  (L)  is  computed  using 
the  equation: 


This  formula  is  the  usual  two-term  servo  lag  approx¬ 
imation  with  acceleration  being  estimated  using  two 
velocity  measurements. 

The  computation  of  regeneration  rates  for  regenera¬ 
tive  sensor  tracking  is  based  on  data  acquired  by  the 
sensor.  This  is  again  in  agreement  with  the  informa¬ 
tion  flow  in  the  real  world  situation.  Estimates  of 
current  rates  (P„)  are  computed  by  extrapolating 
smoothed  rates  computed  at  the  last  data  acquisition 
event.  Thus: 

p;  =  CTOP  [X.A(  +  X.At  At,  *  Y  At  +  Y.At  At,  Z.At 
+-2-At  At.  X.^t’i-Af^-Atl  (3’&) 


The  current  estimate  is  combined  with  the  estimate 
computed  on  the  previous  entry  in  manner  similar  to 
the  computation  for  lag  Thus: 


(3.7) 


A  correlated  sequence  of  random  sensor  error 
(noise)  is  generated  as  follows:  A  three-component 
vector  (N0)  is  computed  using  the  relationship: 


N0=aN.At+dX 


(3.8) 


Where:  X  is  a  three-component  vector  of  inde¬ 

pendent  normally-distributed  variables 
with  mean  zero  and  variance  one. 

k  =  e'A,/-20  (3.9) 

a  =  U-k2)  (3.10) 
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0  =  (l-k2)k  (3.11) 

This  vector  is  then  transformed  into  random  error  in 
range  (R).  azimuth  (9).  and  elevation  (^)  using: 

Nr=N,or  (3.12) 

Ng  =  Nj  o^/(R  cos  d)  (3.13) 

N0-N3ag/R.  (3.14) 


Several  computer  runs  were  made  using  a  more  com¬ 
plex  error  equation.  This  equation  is  presented  in 
Appendix  A. 

The  full  sequence  of  events  in  equation  form  is: 

=  CTOP  irol  +T0  +  Kg0  ♦  N0  (3. 1 5) 


a.  Prediction. 

b.  Gun  Motion  and  Firing. 

c.  Regenerative  Tracking  Switch. 

The  Prediction  and  Gun  Motion  Events  will  he 
scheduled  whenever  the  aircraft  is  visible  to  the  sensor 
and  they  are  not  currently  scheduled.  They  will  be 
unscheduled  whenever  the  aircraft  is  not  visible  and 
they  are  currently  scheduled.  Once  they  have  been 
scheduled,  they  will  continue  to  reschedule  themselves 
at  their  specified  period  until  unscheduled  by  action  of 
the  Sensor  Data  Acquisition  Event. 

The  Regeneration  Switch  ‘ON’  event  will  be  sched¬ 
uled  when  the  aircraft  is  first  sensed  and  the  Regenera¬ 
tive  Switch  ‘OFF’  event  will  be  invoked  whenever  the 
aircraft  becomes  invisible  to  the  sensor  (in  this  case  the 
event  is  not  scheduled  but  invoked  directly  from  the 
Sensor  Data  Acquisition  Event). 

3.2  REGENERATIVE  TRACKING  SWITCH  EVENT 


sensed  values  (3.16) 


(3.17) 


l  smoothed  (3.18) 
I  valuer. 


(3.19) 

On  any  given  run,  the  inclusions  of  lag.  regenera¬ 
tion.  and  noise  are  optional.  Further,  the  number  of 
points  to  be  used  in  smoothing,  as  well  as  the  particu¬ 
lar  coefficients  to  be  used,  are  also  specified  at  input 
time 

Some  further  details  are  worthy  of  note.  First,  since 
regeneration  depends  on  the  presence  of  smoothed 
data  values  for  at  least  two  time  periods,  problems  can 
be  encountered  if  attempts  are  made  to  start  regenerat¬ 
ing  rates  on  the  first  set  of  sensed  data.  Further,  since 
no  provision  is  currently  provided  for  start-up  proce¬ 
dures  in  the  data  smoother,  initial  values  for  smoothed 
data  may  be  unstable  and  result  in  unrealistically  large 
prediction  errors.  (In  a  real  system  the  problems  would 
be  avoided  by  providing  a  start-up  procedure  in  the 
data  smoother )  These  problems  are  avoided  by  delay¬ 
ing  the  first  regeneration  of  rales  and  the  first  predic¬ 
tion  event  by  one  smoothing  time 

The  preceding  comments  lead  to  another  important 
function  of  this  event  routine  This  event  is  responsible 
for  the  initial  scheduling  and  cancelling  of  the  follow¬ 
ing  events: 


C;=PT0C[P^1 
‘j  C-jAt 


n*l  _ 
j*0 


n-1 

-o  ”53  bj  C-jdkt 
j=0 

.  n-1 

£o  =  £  cj  C-j^t 
1=0 


The  function  of  this  event  is  to  turn  the  regeneration 
switch  ‘ON*  and  ‘OFF’  (see  comments  in  the  descrip¬ 
tion  of  the  Sensor  Data  Acquisition  Event).  If  regener¬ 
ation  has  been  specified  in  the  input  and  it  is  invoked 
with  an  ‘ON’  signal,  it  turns  the  switch  ‘ON.’  If  it  is 
invoked  with  an  ‘OFF’  signal  it  turns  the  switch  OFF.’ 
If  regeneration  has  not  been  specified  in  the  input,  the 
invocation  of  this  event  has  no  effect  on  the  simula¬ 
tion. 

3.3  PREDICTION  EVENT 

The  Prediction  Event  models  that  portion  of  the  real 
gun  system  in  which  sensor  information  is  extrapo¬ 
lated  into  the  future  in  order  to  compute  lead  angles 
for  gun  pointing.  This  event  routine  can  operate  in  two 
modes.  The  first,  standard  mode,  makes  no  restraining 
assumptions  on  the  aircraft  trajectory.  The  second, 
defense  of  a  known  point,  assumes  that  a  known  point 
is  the  object  of  attack  and  biases  predictions  according 
to  that  assumption. 

The  sequence  of  computations  which  occur  when  an 
entry  is  made  to  this  event  routine  are  as  follows:  First, 
the  Time  of  Flight  Module  is  invoked  which,  with  the 
aid  of  the  Prediction  Algorithm  Module,  computes  a 
time  of  flight.  The  returned  value  for  time  of  flight  is 
combined  with  current  time  and  the  Prediction  Algo¬ 
rithm  Module  is  invoked  directly  to  obtain  the  re¬ 
quired  predicted  position.  At  this  point,  if  standard 
prediction  has  been  specified,  a  branch  is  executed  to 
the  area  which  stores  data  in  the  PREDICTED  POSI¬ 
TION  File.  If  defense  of  a  known  point  has  been 
specified,  the  process  continues. 

As  mentioned  above,  when  defense  of  a  known 
point  is  specified,  predictions  are  biased.  The  bias  may 
take  two  forms  based  on  either  the  assumption  that  the 
aircraft  will  turn  toward  the  defended  point,  or  the 
assumption  that  the  aircraft  will  follow  a  straight  line 
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from  its  present  position  to  the  defended  point.  Con¬ 
sideration  of  the  geometry  of  the  situation  reveals  that 
neither  of  these  biases  are  appropriate  under  certain 
conditions.  For  this  reason,  even  when  bias  is  specified, 
a  particular  prediction  may  or  may  not  be  biased.  The 
determining  factor  is  the  position  and  velocity  of  the 
aircraft  relative  to  the  defended  point.  If  the  defended 
point  is  in  front  of  the  aircraft,  a  bias  is  considered,  in 
front  of’  should  be  understood  as  an  angle  of  less  than 
90  degrees  between  the  aircraft  horizontal  velocity 
vector  and  the  horizontal  projection  of  the  vector  f rom 
the  aircraft  to  the  defended  point.  Since  the  dot  pro¬ 
duct  of  two  vectors  is  positive,  if  and  only  if  the 
absolute  value  of  the  angle  between  the  vector  involved 
is  less  than  90  degrees,  the  decision  rule  becomes: 

a.  Consider  bias  if:  V.R  >  0. 

b.  Do  not  consider  bias  if :  F.R  £  0. 

Figure  3-1  displays  the  geometry  of  the  decision 
rule.  If  the  defended  point  lies  within  the  shaded  area 
bias  is  considered. 

Whenever  the  defended  point  lies  within  the  bias 
consideration  area  of  the  aircraft,  the  geometry  is 
further  examined  to  determine  if  a  bias  can  be  used, 
and  if  so,  which  of  the  two  biases  is  appropriate.  Three 
situations  are  possible;  (1)  the  aircraft  heading  is  such 
that  a  projection  of  its  current  velocity  vector  would 
pass  within  a  specified  (during  initialization)  distance 
(Rc)  of  the  defended  point,  (2)  a  standard  turn  of 
specified  (during  initialization)  acceleration  followed 
by  a  straight  line  flight  path  would  cause  the  aircraft  to 
pass  over  the  target,  or  (3)  neither  of  the  above. 

The  first  situation,  namely  the  straight  line  bias,  is 
the  easiest  to  explain  (see  Figure  3-2).  The  heading 
angle  change  0  is  computed  as: 


The  closest  approach  to  the  defended  point  is  ap¬ 
proximated  by  the  arc: 


00634  301 


Figure  3-1  Bias  Consideration  Region 


POINT 
COP  T4  XI? 

Figure  3-2  Straight-Line  Bias  Situation 

S=R0  (3.21) 

Linear  bias  is  used  if  S  s  R,. 

The  linear-bias  computations  are  accomplished  in 
the  Time  of  Flight  Module  and  the  Prediction  Algo¬ 
rithm  Module.  If  a  situation  is  such  that  linear  bias  is 
appropriate,  an  indicator  switch  is  set  and  point  to 
point  interpolation  will  be  used  in  the  Prediction 
Algorithm  Module. 

If  the  situation  is  such  that  linear  bias  is  not  appro¬ 
priate,  a  check  is  made  to  see  if  a  turn  followed  by  a 
straight-line  path  can  bring  the  aircraft  over  the  target. 
This  check  is  made  as  follows:  The  radius  of  turn  is 
computed  using: 

6  =  V2/ng  (3.22) 

where: 

V  =  the  aircraft  velocity 

ng=  the  standard  turn  acceleration  (n  is  input) 

The  heading  and  turn  radius  are  used  to  find  the 
center  of  curvature  of  the  turn  and  the  distance  (L) 
from  this  point  to  the  defended  point.  Whenever  the 
turn  radius  is  less  than  the  distance  from  the  center  of 
curvature,  a  turn  is  possible  (Figure  3-3).  If  no  turn  ts 
possible,  we  have  situation  three:  that  is,  no  bias  can 
be  used  (Figure  3-4).  Under  such  circumstances  a 
branch  to  the  data  storage  area  is  executed  and  stan¬ 
dard  values  for  prediction  are  used. 

If  a  turn  followed  by  a  linear  segment  is  to  be  used, 
the  process  continues.  Figure  3-5  presents  the  d'  iled 
geometry  of  the  turn.  The  new  heading  (0)  anc  the 
time  to  turn  (Tt)  are  computed  using  the  following 
formulas: 

Rp  =  v/R02  -  26  R0sinH  (3.23) 
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Figure  3-3  Turn  to  Attack  Bias  Situation  <S  <  L) 


(R0  -  5  sin  H)  -  R.  5  cos  H 


\»Xq  u  Jill  »•/  ’  1-Vp  o  VVJ  II 

j  (3.24) 

Rp  (RQ  -  6  sin  H)  +  62cosH 

Afl  =  A  +  H 

(3,25) 

'■If 

(3.261 

tt  =  e/6 

(3.27) 

e  -  e0  +  &e 

(3,28) 

eluded  The  time  of  flight  to  (he  end  point  of  the  turn 
is  compared  to  the  time  to  turn.  If  the  time  of  flight 
exceeds  the  time  to  turn,  the  biased  prediction  is  used, 
otherwise  a  branch  to  the  storage  area  is  made  and  a 
standard  prediction  results. 

When  turn  bias  is  included,  values  for  the  velocity  in 
the  x  and  y  direction  are  computed.  The  prediction  in 
these  directions  is  a  linear  extrapolation  along  the  line 
from  the  end  of  the  turn  to  the  defended  point.  In  this 
case,  the  Z  coordinate  is  unbiased. 

For  all  cases,  the  last  action  is  a  branch  to  the  data 


Figure  3-5.  Geometry  of  a  Turn  Followed  by  a 
Linear  Segment 


storage  area  where  the  predicted  position  is  trans¬ 
formed  to  polar  coordinates  and  stored  along  with  the 
estimated  time  of  closest  approach  as  the  most  recent 
entry  in  the  PREDICTED  POSITION  File. 

Figure  3-6  summarizes  the  logical  structure  of  the 
Prediction  Event 

3.4  GUN  MOTION  AND  FIRING  EVENT 

The  function  of  this  event  is  to  compute  the  closest 
approach  of  a  shell  fired  at  current  clock  time  given: 
the  history  of  predicted  positions,  the  ballistic  coeffi¬ 
cients,  and  the  actual  aircraft  flight  path.  Gun  servo 
lag  and/or  servo  regeneration  can  be  included.  When 
included,  lag  and/or  regeneration  characteristics  are 
determined  by  input  parameter  specification. 

In  the  following  discussion,  upper  case  letters  are 
used  to  designate  predicted  values,  subscripts  designate 
the  time  to  which  the  values  pertain,  primed  upper 
case  letters  are  used  to  indicate  actual  gun  pointing 
angles,  and  the  letters  L  and  R  are  used  to  designate 
Lag  and  Regeneration  respectively.  Standard  dot  nota¬ 
tion  is  used  to  indicate  rates  and  accelerations.  In  this 
notation,  the  azimuth  and  elevation  become 

*o  =  eo  +  Lo  +  Ro  <3-29> 


m 


i 


i 


\ 


1 

2 


3-4 


*c-*o+Lo  +  Ro 


(3.30) 


COMPUTE  \EY, 
P«EOiCT>ON 
USING 
fNE  8«AS 


The  lag  computation  is  based  on  the  four  latest 
predicted  positions.  These  four  values  are  used  to 
estimate  current  rate  and  acceleration  values  for  the 
coordinate  in  question,  which  is  in  turn  used  in  the 
expression: 


L 


A 


Kv  Kg 


(3.31) 


When  the  estimates  are  substituted  into  the  above 
expression  the  lag  equations  become: 


1.0 

2.0 

|  Ao  + 

1-1.0  -5.0 

KvAt  + 

M’2 

|  K2At2 

4.0 

-!.0 

v,2i 

A-2At  +  | 

KAt2 

'-3At 


'-At 

(3.32) 

(3.33) 


Values  for  regeneration  are  computed  as  follows: 
The  current  predicted  position  is  transformed  to  carte¬ 
sian  coordinates.  These  are  combined  with  current 
values  for  smoothed  cartesian  rates  and  transformed 
back  to  polar  coordinates.  The  resulting  polar  rates  are 
used  to  estimate  polar  acceleration  using  the  following 
equations: 


svv+V 

(3.34) 

S  =  (XpXp  +  YpYp)/S 

(3.35) 

it  =  -:.os<i/D 

(3.36) 

-2.0  D  DO  -  S0Z_ 

.  •  P 

(3.37) 

These  values  are  used  in  the  standard  two-term  lag 
estimate  yielding. 


R0  = 


0  |  l 

•V*1  KaAt2 


(3.38) 
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the  Prediction 


V 


(3.39) 


The  projectile  position,  at  current  time  plus  time  of 
flight  (estimated  time  of  closest  approach),  is  given  by: 
the  predicted  range  <D,i.  the  adjusted  azimuth  id \). 
and  the  adjusted  elevation  «f>,)  The  projectile  velocity 
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tc  =  |RxV|/V- V 


(3.41) 


is  assumed  to  he  in  the  radial  direction  with  magnitude 
given  by. 

.1 

vp(t) "  (nr  vma'' 1 + 1-°t 177  vm  <3-40> 

where: 

l  =  time  of  flight. 

Vm  »  muzzle  velocity. 

a,  K  =  ballistic  coefficients  from  the  Time  of  Flight 
module. 

Projectile  position  and  velocity  are  transformed  into 
cartesian  coordinates  for  the  computation  of  the  actual 
closest  approach  values. 

The  time  of  closest  approach  is  computed  by  correct¬ 
ing  the  estimated  time  of  closest  approach  which  is 
computed  by  the  predictor  The  correction  is  made  by 
assuming:  that  the  aircraft  and  projectile  are  moving 
at  constant  rates  for  the  time  period  necessary  for  the 
projectile  to  get  from  either  the  predicted  position  to 
position  of  closest  approach,  or  in  the  case  when  the 
closest  approach  occurs  earlier  than  the  predicted  time, 
the  time  necessary  to  travel  from  the  closest  approach 
to  the  predicted  position.  Once  this  assumption  is 
made,  the  time  from  predicted  position  to  closest 
approach  can  be  computed  as: 


where: 

R  =  position  of  aircraft  relative  to  the  projectile 
at  predicted  impact. 

V  =  velocity  of  the  aircraf  t  relative  to  the  projectile 
at  the  predicted  impact. 

The  time  of  closest  approach  is  then  given  by: 

t  =  tp+tc  (3.42) 

The  position  of  the  shell  at  closest  approach  is  given 
by: 

WVslc  (343) 

If  the  error  components  are  desired,  they  are  com¬ 

puted  as: 

E=R  +  Vtc  (3.44) 

With  R,  V,  and  tc  defined  as  above. 

The  time  of  fire,  the  range,  and  the  miss  distance  at 
closest  approach  are  stored  in  the  SHELL  POSITION 
File.  There  are  several  versions  of  this  event  routine 
which  differ  in  the  contents  stored  in  the  file  for 
printing. 
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SECTION  4 
UTILITY  PACKAGE 


4.1  AIRCRAFT  POSITION  ROUTINE 

The  Aircraft  Position  Routine  provides  the  position 
and  velocity  of  the  aircraft  in  cartesian  coordinates 
centered  at  the  gun  for  the  simulation  time  specified  in 
the  invoking  command 

Paths  are  represented  as  sequences  of  continuous 
path  segments  proceeding  from  an  initial  point  of 
known  position  and  velocity.  Each  path  segment  is 
either  a  segment  of  a  line,  a  segment  of  a  circle,  or  a 
segment  of  a  logarithmic  spiral.  All  path  segments  are 
in  the  horizontal  plane  coupled  with  motion  of  con¬ 
stant  acceleration  (possibly  zero)  in  the  vertical  direc¬ 
tion.  Smooth  positions  and  velocities  are  assured  by 
using  the  end  conditions  of  one  segment  as  the  initial 
conditions  of  its  successor.  The  end  point  of  a  segment 
is  determined  relative  to  its  initial  point  with  the 
‘duration’  of  the  segment  being  specified  by:  length  for 
a  line  segment,  and  degrees  of  turn  for  both  circular 
and  spiral  segments.  The  only  position  and  velocity 
specification  made  is  for  the  initial  point.  Further  path 
description  is  in  the  form  of  the  following  statements: 

a  Line  ‘X‘  meters  in  length,  tangential  acceleration 
'At,'  and  vertical  acceleration  ‘A2.’ 

b.  Circular  turn  '4>'  degrees  right  or  left  radial 
acceleration  ‘A«’  and  vertical  acceleration  'Ai.' 

c.  Spiral  turn  ‘<f>’  degrees  right  or  left,  radial  accel¬ 
eration  ‘An,’  tangential  acceleration  ‘At,’  and 
vertical  acceleration  ‘A?.’ 


For  linear  segments  position  and  velocity  are  deter¬ 
mined  according  to  the  equations: 


XfO’Xo-^Vot  +  ^ATt2^cos0 

(4,1) 

Y(t)  =  Y0t(V0t  +  iATt2)sin0 

(4.2) 

Z(t)  =  Z0  +  VZot  ♦  jAZt2 

(4.3) 

Vx(t)  =  -(V0  +  Ajt)  cos0 

(4.4) 

Vz(0  =  VZo*AZt  (4.6) 

Where; 


0  “  aircraft  heading 

For  circular  segments,  position  and  velocity  are 
determined  according  to  the  equations: 


R  =  V02/AR  (4.7) 

6  =  V0/R  (4.8) 

0(t)  =  CQ  +  St  (4.9) 

+  R  [cos  0O  -  cos  0(i))  (4.10) 

Y(t)=Y0  +  R  [sin  0O  +  sin0(t)|  (4.11) 

Z(t)  =  Z0  +  VZot  +  1azi2  (4.12) 

Vx(t)=  -V  sin  0(t)  (4.13) 

Vyft)  =  V  COS  0(t)  (4.14) 

VZt‘)  =  VZo  +  AZt  (4.15) 

For  spiral  segments  position  and  velocity  are  deter¬ 
mined  according  to  the  equations: 

X=At/Ar  (4.16) 

V(t)*V0  +  Ajt  (4.17) 

R  =  Vo2/Ar  (4.18) 


0(t)  =  0o  -  Rlogg  [V2(t)/2XV02] 


xct)  =  X0  -  »  j  2x|sin  0O  -  ^  sin  <0(t)) 

Ar  0+4X2)  I  L  V02  J 

-  j^cos0o-  ^11  cos (0(1))]  |  (4.20) 

Y(t)  =  Y0  =  -— ^  -  hxf^  cos0(,)-cos0ol 

Ar  ( I  +  4X2)  |  (_v02 


COS  0(t)  -  cos  0. 


Vy(t)  =  (vo  ♦  Aft)  sin  0 


(4.21) 


(4.22) 


«»*Z0  +VZot  ♦  JaZi2 

V x(t)  *  ~V(t)  sin  10<U| 

Vy(t)  =  V(t)  cos  10(01 

VZ«>  =  V4  +  Azt 


(4.23) 

(4.24) 

(4.25) 


‘start  ■'a1 


Ra(0) 


fca(0) 


v  0)  +  Vm2-Va2(0)  /(vm2-va2(0)) 


(4.28) 


where: 


4.2  TIME  OF  FLIGHT  MODULE 

The  purpose  of  the  Time  of  Flight  Module  is  to 
compute  the  time  required  for  a  projectile  to  travel 
from  the  gun  to  a  future  aircraft  position.  The  prob¬ 
lem  is  somewhat  complicated  by  the  fact  that  the 
predicted  position  is  dependent  on  the  time  which  is  of 
course  based  on  the  time  of  flight.  Closed-form  solu¬ 
tions  of  the  problem  are  possible  only  under  some 
rather  restrictive  assumptions.  For  this  reason  numeri¬ 
cal  methods  have  been  employed.  The  time  of  flight 
computation  is  based  on  the  fact;  that  at  the  time  of 
predicted  impact,  the  range  to  the  aircraft  (R,)  and  to 
the  projectile  (R*)  must  be  equal.  In  essence  then,  the 
Time  of  Flight  Module  numerically  solves  the 
equation: 

f(t)  =  Ra(t)  -  Rp(t)  =  0  (4.26) 

for  the  smallest  value  of  t. 

The  values  of  aircraft  range  are  provided  by  the 
Prediction  Algorithm  Module  and  are  dependent  on 
the  values  currently  stored  in  the  SENSED  DATA  File 
as  well  as  the  particular  options  for  prediction  being 
simulated.  Values  for  R„(t)  are  computed  using  the 
equation: 

Rp(t)  *  a[(bt  +  1  )c  —  1 1  (4.27) 

where: 

a  =  1.0  K(cl-  2)  Vma-; 

b  =  K(a  -  1 )  Vma' 1 

c  =  (a  -  2)/(a  -  1) 

a  =  input  variable 

K  =  input  variable 
Vm  =  input  variable 

A  start  value  for  the  solution  is  obtained  by  assum¬ 
ing  a  projectile  which  does  not  decelerate,  and  an 
aircraft  which  is  in  unaccelerated  motion: 


Vm  =  muzzle  velocity  of  gun. 

Va  =  current  aircraft  velocity. 

ka  *  current  range  rate  of  the  aircraft. 

An  initial  computation  of  the  function  f(t)  *=  R,(t)  - 
R,(t)  is  made  using  the  start  value.  The  process  is 
continued  with  each  successive  computation  being 
made  at  a  new  time  which  is  the  current  value  incre¬ 
mented  by  a  function  of  the  start  value. 

Once  a  sign  change  has  been  obtained,  a  combina¬ 
tion  of  the  secant  method  and  interval  bisection  is 
employed  until  the  desired  accuracy  is  obtained.  In 
general,  no  more  than  five  compulations  are  required 
to  derive  the  required  accuracy. 

4.3  PREDICTION  ALGORITHM  MODULE 

The  function  of  the  Prediction  Algorithm  Module  is 
to  extrapolate  current  available  aircraft  position  and 
velocity  data  to  a  specified  future  time  point;  thereby, 
producing  an  estimate  of  aircraft  position  for  the 
specified  time.  The  module  has  the  capability  of  per¬ 
forming  the  extrapolation  by  a  number  of  methods; 
each  of  which  is  somewhat  flexible  in  that  its  specific 
behavior  is  determined  by  method  parameters.  The 
specific  method,  employed  during  a  simulation  run.  is 
determined  by  the  designation  of  the  general  method 
type  and  its  required  parameters  during  the  initializa¬ 
tion  phase  of  the  run.  The  position  data  on  which  the 
module  operates  are  all  derived  from  the  output  of  the 
sensor  module.  These  data  are  accessed  via  the  commu¬ 
nications  pool  (specifically  the  SENSED  DATA  File). 

Current  module  capabilities  include  three  general 
methods: 

a.  Quadratic  (includes  linear  by  parameter  specifi¬ 
cation). 

b.  Simple  Polar. 

c.  Defense  of  a  Known  Point  (in  combination  with 
Quadratic). 

4.3.1  Quadratic  Prediction 

The  quadratic  prediction  is  made  using  smoothed 
position  (X.  Y.  Z).  velocity  (X.  Y,  Z).  and  acceleration 
(X,  Y.  Z)  values  computed  and  stored  in  the  SENSED 
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POSITION  File  by  the  TRACK  Module.  The  predic- 
tion  is  computed  as:  a  linear  estimate  corrected  for 
accelerations  in  the  horizontal  plane,  and  changes  in 
total  aircraft  speed  induced  by  changes  in  aircraft 
altitude.  The  basic  estimate  for  a  time  increment  r  is 
given  by: 


X  ■  «  X  +  Xr 

(4.29) 

Y  -  Y  +  Vt 

(4.30) 

Z  =  Z  +  It 

(4.31) 

The  correction  for  acceleration  in  the  horizontal 
plane  is  made  whenever  the  absolute  value  of  the 
estimated  heading  rate  6  of  the  aircraft  exceeds  the 
value  of  an  input  threshold  parameter.  The  specifica¬ 
tion  of  an  excessively  large  threshold  value  leads  to 
suppression  of  the  acceleration  correction  and  linear 
prediction  results.  The  heading  rate  is  estimated  by  the 
equation: 

i  =  (xi  -  XY )/(X2  +  Y2)  (4.32) 


When  the  correction  is  included,  the  resultant  pre¬ 
dictions  for  X  and  Y  are: 


X  =  X  +  Xf  - 


Y '  *  Y  +  Xt  + 


‘i(r2  +  rAt+^; 

)v 

(4.33) 

U(r2  +  TAt  +  ^r, 

)x 

(4.34) 

The  correction  for  change  in  aircraft  speed,  due  to  a 
change  in  altitude,  is  made  whenever  the  absolute 
value  of  the  aircraft  velocity  in  the  vertical  direction 
exceeds  the  value  of  an  input  threshold  parameter.  As 
in  the  case  of  the  acceleration,  the  correction  can  be 
suppressed  by  specification  of  an  excessive  value  for 
the  parameter.  When  the  correction  is  included,  the 
resultant  predictions  are: 

!LY+ g(XZ/v2)) 

(4.35) 

[0X-g(YZ/v2)) 

(4.36) 

[gfZ2/V2)j 

(4.37) 


The  Simple  Polar  Prediction  is  made  by  estimating 
the  range  rate  R.  the  azimuth  rate  A.  the  elevation  rate 


X'  =  X  +  Xr  -  j^  +  TAt-t-^) 

Y '  =  Y  +  Yt  +  ^(r2  +  tA  t  + 

z '  *  z  +  Zt  +  +  +  TT") 

where: 

At  =  Sample  Interval 
V  =  Estimated  Aircraft  Velocity 
g  =  Acceleration  due  to  gravity 

4.3.2  Simple  Polar  Prediction 


E;  assuming  constant  rates,  and  using  linear  extrapola¬ 
tion.  Polar  coordinate  rates  are  estimated  by  trans¬ 
forming  the  cartesian  coordinates  of  the  last  two  en¬ 
tries  in  the  SENSED  POSITION  File  and  computing 
first-divided  differences.  Thus: 


R  H  (Rt  ~  Bt-AtWAt 

(4.38) 

A  =  (At  -  A ,.£,)/ At 

(4.39) 

1  «  (Et  ~  Et-At )M‘ 

(4.40) 

Estimated  range,  azimuth,  and 
by: 

elevation  are  given 

R  *  Rt  +  Rt 

(4.41) 

A  '  *  At  +  At 

(4.42) 

E  °  Et  +Et 

(4.43) 

The  required  cartesian  estimates  are  obtained  by 
transforming  R’.  A’,  and  E’. 

4.3.3  Defenae  of  a  Known  Point 


The  majority  of  the  logic  and  computation  for  the 
defense  of  a  known  point  option  are  in  the  Prediction 
Event  Module.  However,  for  time  of  flight  computa¬ 
tions,  use  is  made  of  the  Prediction  Algorithm  Module. 
During  such  entries  to  the  module,  one  of  two  modes 
of  prediction  are  possible:  Quadratic  (discussed  previ¬ 
ously)  or  point-to-point  linear  interpolation.  An  indi¬ 
cator  switch  set  in  the  Prediction  Event  Module  con¬ 
trols  which  method  is  used  on  any  given  entry  to  the 
module.  The  following  pertains  to  point-to-point  linear 
interpolation  used  for  a  portion  of  the  attack  of  a 
known  point  option. 

Figure  4-1  presents  the  geometry  for  the  point-to- 
point  linear  interpolation.  The  point  Q  represents  the 
last  sensed  aircraft  position  and  P  the  point  being 
defended.  The  prediction  method  assumes  that  the 
aircraft  will  travel  along  the  line  connecting  the  points 
at  rates  which  are  consistent  in  magnitude  to  the 
current  smoothed  aircraft  speed.  Distance  along  the 
line  is  computed  using  the  equations: 


AX  =  Xp  -  Xq 

(4.44) 

AY  =  Yp  -  Yq 

(4  45) 

AZ  =  Zp  -  Zq 

(4.46) 

D  =  \J  AX2  +  AY-  +  AZ2 

(4.47) 
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(4.52) 


Ut)  * 

Vt  + 

;g(t2 

+  tTs  + 

Ts2) 
6  / 

'  D 

(4.48) 

The  desired 

estimate  is 

then 

obta 

ined  by 

interpolation 

X  = 

Xq 

+  AX 

Ut).D 

(4.49) 

Y  = 

yQ 

+  AY 

■  L(t)/D 

(4.50) 

Z  = 

Zq 

+  AZ 

■  Ut)/D 

(4.51) 

4.4  RANDOM  NUMBER  GENERATOR 

The  function  of  the  Random  Number  Generator  is 
to  produce  samples  of  random  variables  from  either  a 
population  uniformly  distributed  on  the  interval  [0,1], 
or  from  a  population  which  is  normally  distributed 
with  mean  zero  and  variance  one.  The  uniform  ran¬ 
dom  variables  are  generated  by  using  an  additive 
congruous  process  with  31  seed's.  This  method  was 
based  on  work  which  was  done  at  MIT’s  Lincoln 
Laboratory  in  1958  by  B.F.  Green,  J.E.  Smith,  and 
L.  Klem.  The  sequence  of  samples  is  produced  by 
using  the  following  equations- 


l 
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Figure  4-1  Geometric  Representation  of  the  Point 
to-Point  Linear  Interpolation 


Xj  »  (Xj.i  +  Xj.n)  Mod  2r  j>n 

Rj-n+1  ■  Xj/2r  (4.53) 

where : 

n  =  31 
r«  29 

Xj  =  an  input  seed  variable 
X2  -  X31  ■  prestored  seed  values 

Rj  ■  the  required  uniformly  distributed  random 
numbers. 

Normally-distributed  variables  are  computed  from 
the  uniformly-distributed  random  variables  using  the 
wel!  known  Central  Limit  Theorem.  The  exact  equa¬ 
tion  employed  is: 

12 

N  R,  -  6.0  (4.54) 

i  =1 

where; 

N  =  the  required  sample 

Rj  =  uniformly  distributed  random  variables 

This  module  has  the  ability  to  return  any  number  of 
normal  or  uniform  numbers  on  any  given  entry.  The 
sequence  of  samples  generated  is  a  function  of  the  first 
seed  value  which  is  an  input  variable. 

4.5  COORDINATE  TRANSFORMATION 

There  are  two  coordinate  transformation  programs. 
One  converts  polar  coordinates  to  cartesian  coordinates 
and  the  other  converts  cartesian  coordinates  to  polar 
coordinates.  Both  programs  convert  the  three  position 
coordinates  plus  their  respective  rates.  The  transforma¬ 
tions  are  not  standard  due  to  the  fact  that  the  azimuth 
angle  is  measured  counter-clockwise  from  the  positive 
y-axis.  This  system  was  used  in  an  attempt  to  remain 
as  compatible  as  possible  with  systems  employed  in  the 
University  of  Michigan  Gun  Model. 

The  cartesian  to  polar  transformation  uses  Equations 
4.55  through  4.61.  The  equations  are  computed  in  the 
order  shown  and  earlier  results  are  used  in  later  com¬ 
putation  in  order  to  minimize  computer  time. 

R  =  >/x2  +  Y2  +  Z2  (4.55) 

S  =  \/x2  +  Y2  (4.56) 
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0  =  Tan->  (-X/Y) 

(4.57) 

<t>  =  Tan'1  (Z/S) 

(4.58) 

R  =  (XX  +  YY  +  ZZ)/R 

(4.59) 

0  =  (XY  -  XY)/S2 

(4.60) 

0  =  (Z  -  ZR/R)/S 

(4.61) 

The  polar  to  cartesian  conversion 
equations: 

uses  the  following 

S  ■  R  cos  <t> 

(4.62) 

Z*  Rsin  <j> 

(4.63) 

X  ■  -S  sin  8 

(4.64) 

Y  *  S  cos  0 

(4.65) 

t»R/r-z£/s 

(4.66) 

*-XT-Ytf 

(4.67) 

t-Xd  +  YT 

(4.68) 

t  -  Zft/R  +  S<fr 

(4.69) 

4.6  REPORT  GENERATORS 

Report  generators  are  provided  to  compute  and 
output  measures  or  effectiveness  for  the  gun  system 
being  simulated.  Due  to  the  variability  of  the  measures 
required  for  different  aspects  of  the  system  study,  two 
general  types  of  report  generators  have  been  provided. 
For  detailed  analyses  of  the  effect  of  a  single  aspect  of 
the  system,  periodic  reports  are  provided.  For  the 
comparison  of  the  effectiveness  of  the  system  as  a 
whole,  summary  reports  are  produced  by  an  exit  report 
generator. 

4.6.1  Periodic  Reports 

Periodic  Reports  are  generated  by  special  report 
events.  They  may  be  either  scheduled  or  invoked  by 
another  module.  In  general,  they  have  access  to  any 
data  in  the  communication  pool  but  may  not  alter 
these  values.  Figure  4-2  presents  an  example  of  the 
output  from  a  periodic  report. 

This  report  presents  a  description  of  the  aircraft  and 
shell  positions  at  the  time  of  closest  approach  of  the 
shell  to  the  aircraft.  The  time  entries  include:  the  clock 
time  when  the  shell  was  fired,  the  estimated  time  of 
flight,  and  the  clock  time  at  the  closest  approach. 


The  first  line  presents  shell  range  in  meters,  azimuth 
in  mils,  and  elevation  in  mils.  The  second  line  presents 
aircraft  range  in  meters,  azimuth  in  mils,  and  elevation 
in  mils.  The  third  line  presents  the  difference  in  posi¬ 
tion  (measured  in  meters)  in  the  range,  azimuth  and 
elevation  directions.  The  fourth  line  presents  the  errors 
in  the  azimuth  and  elevation  coordinates  measured  in 
mils.  This  particular  report  was  generated  at  a  time 
interval  of  0.20  second. 

4.6.2  Exit  Report  Generator 

As  stated  previously,  the  exit  report  generator  is 
concerned  with  the  performance  of  the  system  as  a 
whole.  The  measures  of  effectiveness  which  it  produces 
are  therefore  biased  towards  integrations  of  perform¬ 
ance  over  the  entire  duration  of  the  course.  Other 
measures  may  be  produced  and  outputted  by  indicat¬ 
ing  the  requirement  to  do  so  on  a  print  options  card. 

With  no  output  indicated,  the  exit  report  will  gener¬ 
ate  for  each  replication: 


o 


where: 


pss  =  single  shot  kill  probability 
T  =  duration  of  the  aircraft  flight 

In  addition,  at  the  end  of  the  last  replication  the 
mean  and  standard  deviation  of  K,  over  the  replication 
is  output.  Furthermore,  the  distribution  over  the  repli¬ 
cates  of  the  average  pss  over  bursts,  and  the  distribu¬ 
tion  of  miss  distances  for  single  samples  are  output. 
Figure  4-3  presents  an  example  of  the  summary  pro¬ 
duced  by  the  exit  report. 

Information  may  also  be  output  at  the  individual 
replicate,  burst,  or  sample  level.  Figures  4-4  through 
4-7  present  examples  of  the  outputs  at  these  various 
levels. 

Additional  Outputs 

In  any  particular  study  for  which  the  current  outputs 
are  not  suitable,  additional  reports  can  be  generated 
with  minor  modification  to  the  model.  In  most  cases, 
the  modifications  will  be  the  replacement  of  an  existing 
report  subroutine  with  one  more  appropriate  to  the 
study  being  conducted. 
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SECTION  i 

COMMUNICATION  POOL 


The  parameter  values  stored  in  the  communication 
pool  can  be  viewed  as  three  data  files  and  a  switch  hie. 
The  data  files,  which  include  SENSED  POSITION, 
PREDICTED  POSITION,  and  SHELL  POSITION,  are 
used  to  organize  related  data  and  to  restrict  access  to 
those  modules  which  have  a  ‘right’  to  access  the  data  in 
question.  The  switch  file  is  a  collection  of  program 
controlled  parameters  which  are  used  to  alter  the 
logical  flow  of  the  model  and/or  the  system  under 
study. 

5.1  SENSED  POSITION  FILE 

The  content  of  the  SENSED  POSITION  File  is 
presented  in  Table  V-l.  Entries  in  this  file  are  com¬ 
uted  and  stored  by  the  Sensor  Data  Acquisition 
vent.  They  may  be  accessed  only  by  the  Prediction 
Algorithm  Module,  the  Prediction  Event,  or  a  report 
generator. 

6.2  PREDICTED  POSITION  FILE 

The  content  of  the  PREDICTED  POSITION  File  is 
presented  in  Table  V-2.  Entries  in  this  table  are  com¬ 
puted  and  stored  by  the  Prediction  Event.  They  may  be 
accessed  by  the  Gun  Motion  and  Firing  Event  or  a 
report  generator. 

S  3  SHELL  POSITION  FILE 

The  SHELL  POSITION  File  is  filled  by  the  Gun 
Motion  and  Firing  Event  and  accessed  by  the  summary 
report  generator.  The  corner  of  this  file  is  variable 
since  there  exists  several  vi  os  of  the  summary 
report  generator  with  corresponding  modified  versions 
of  the  Gun  Motion  and  Firing  Event.  There  are, 
however,  some  general  principles  that  can  be  stated. 
The  file  will  always  contain  one  entry  for  each  occur¬ 
rence  of  the  Gun  Motion  and  Firing  Event.  The  entries 
will  be  ordered  by  time  of  occurrence.  Entries  for 
which  no  firing  occurred,  due  to  the  aircraft  being  out 
of  range,  will  be  indicated  by  a  negative  sign  on  at 
least  one  of  the  time  parameters  in  the  entry.  Typical 
output  reports  were  presented  in  Section  4.  Examples 
of  entries  in  the  file  which  have  been  implemented  to 
generate  these  reports  include  the  following: 

a.  Time  of  Fire. 

Time  of  Closest  Approach. 

X  of  the  shell  at  Closest  Approach. 

Y  of  the  shell  at  Closest  Approach. 

Z  of  the  shell  at  Closest  Approach. 


b.  Time  of  Fire. 

Time  of  Closest  Approach. 

Error  in  X  at  Closest  Approach. 

Error  in  Y  at  Closest  Approach. 

Error  in  Z  at  Closest  Approach. 

c.  Time  of  Fire. 

Range  at  Closest  Approach. 

Miss  Distance  at  Closest  Approach. 

5.4  SWITCH  FILE 

There  are  five  switches  in  the  SWITCH  File.  Their 
function  is  to  communicate  the  existence  of  situations 
which  require  other  than  normal  or  main-flow  type 
computations.  The  five  switches  and  their  functions 
are: 

a.  ION.  ION  indicates  the  status  of  the  PREDIC¬ 
TION  Event.  If  the  PREDICTION  Event  is  cur¬ 
rently  scheduled,  it  has  value  one.  If  it  is  not 
scheduled.  ION  is  zero.  The  value  of  ION  is  both 
set  and  used  in  the  DATA  ACQUISITION 
Event.  It  is,  however,  set  and  used  on  different 
entries. 

b  REGENTRACK.  This  switch  indicates  whether 
or  not  regeneration  should  be  used  in  the  track¬ 
ing  computation.  Its  value  is  set  by  the  Regenera¬ 
tive  Tracking  Switch  Event  and  used  bv  the 
Sensor  Data  Acquisition  Event. 

c  REGEN-GUN.  This  is  an  ‘on-ofT  switch  to  con¬ 
trol  regenerative  gun  tracking.  Its  value  is  set  at 
system  initialization  time,  and  remains  constant 
throughout  a  simulation  run.  It  is  used  in  the 
Gun  Motion  and  Firing  Module 

d.  DKP-ELAG.  This  is  an  ‘on-olT  switch  for  the 
Defense  of  a  Known  Point  method  of  prediction. 
Its  value  is  set  at  system  initialization  time  and 
remains  constant  throughout  a  simulation  run. 
The  switch  value  determines  the  logical  flow 
within  the  Prediction  Event  Module 

e.  JON.  JON  indicates  whether  or  not  the  point-to- 
point  interpolation  scheme  should  be  used  in  the 
Prediction  Algorithm  Module  Its  value  is  set  by 
the  Prediction  Event  Module  This  switch  will  be 
“olf  unless  linear  bias  is  being  used. 
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APPENDIX  A 

SPECIAL  TRACKER  MODIFICATION 


During  the  course  of  ihe  study,  it  became  apparent 
that  the  tracking  aim-wander  error  equation  should  be 
expanded  to  include  the  capability  to  model  systematic 
errors  and/or  certain  unique  characteristics  of  the 
performance  of  human  trackers.  The  form  of  the 
desired  equation  was  known  but  due  to  a  lack  of  funds 
it  could  not  be  integrated  into  the  model.  A  temporary 
■patch’  was  made  which  enabled  the  model  to  perform 
the  desired  computation  for  a  fixed  set  of  parameter 
values.  A  recompilation  of  the  model  was  required  for 
each  parameter  change. 

The  rull  equations  for  aim  wander  and/or  system¬ 
atic  errors  used  were: 

eD  =  +  Xquu  (A-1) 

Ea  =  UA  +  XAloA  +  X(LagA- RegA)  ■  Dcos(E)] 

(A. 2) 

Ee  =  UE  +  XE  [oE  +  X(LagE  -  RegE)  D|  (A.3) 
Where. 

UD  =  systematic  error  in  range 
UA  =  systematic  error  in  azimuth 
UE  =  systematic  error  in  elevation 


Xq  =  Normal  Random  variable,  mean  “  0, 

variance  =1.  ; 

XA  *  Normal  Random  variable,  mean  *  0, 

variance  =  I . 

XE  =  Normal  Random  variable,  mean  =  0, 
variance  =  I. 

X  =  a  parameter  which  expresses  the  degree  of 
dependence  of  aim  wander  error  on 
coordinate  lag. 

LagA  =  lag  in  azimuth 

Lagg  =  lag  in  elevation 

RegA  =  regenerated  azimuth  rate 

Regg  =  regenerated  elevation  rate. 

Full  integration  of  the  equation  into  the  model 
would  enable  further  study  in  the  area  of  human 
tracking.  It  would  also  open  up  the  possibility  of 
studying  the  effect  of  ihe  transition  from  one  tracking 
mode  to  another  during  an  equipment  malfunction.  A 
failure  could  be  modeled  by  a  special  event  which,  on 
occurrence,  would  modify  the  values  of  appropriate 
parameters  in  the  error  equation.  Thus,  a  change  in 
value  from  zero  to  a  value  appropriate  to  manual 
tracking  for  IV  UA.  L't.  and  X  would  simulate  a 
transition  from  radar  to  manual  tracking. 
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